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r-| ■ ABSTRACT 

\[ We systematically examine the fundamental frequencies of shear torsional oscillations 

in neutron star crusts in a manner that is dependent on the parameter L characterizing 
$H \ the poorly known density dependence of the symmetry energy. The identification of 

c/5 ■ the lowest quasiperiodic oscillation (QPO) among the observed QPOs from giant flares 

in soft-gamma repeaters as the £ = 2 fundamental torsional oscillations enables us to 
constrain the parameter L as L ^ 47.4 MeV, which is the most conservative restriction 
on L obtained in the present work that assumes that the mass and radius of the flaring 
neutron stars range 1.4-1.8 Mq and 10-14 km. Next, we identify one by one a set of the 
low-lying frequencies observed in giant flares as the fundamental torsional oscillations. 
The values of L that can reproduce all the observed frequencies in terms of the torsional 
oscillations coupled with a part of dripped neutrons via entrainment effects are then 
■ constrained as 101.1 MeV ^ L ^ 131.0 MeV. Alternatively, if only the second lowest 

\ frequency observed in SGR 1806—20 has a different origin, one obtains relatively low 

(f) ■ L values ranging 58.0 MeV ^ L ^ 85.3 MeV, which seem more consistent with other 

empirical constraints despite large uncertainties. 
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1 INTRODUCTION 

Neutron stars are believed to form as stellar remnants of core-collapse supernovae, which occur at the end of massive star 
evolution. The density of matter inside a neutron star can easily become higher than normal nuclear density, a situation hard 
to realize in the laboratory. Due to such high densities and resulting strong coupling nature, the equation of state (EOS) for 
neutron star matter remains to be fixed. Consequently, neutron stars act as a unique laboratory to constrai n the properties 



of su ch dense matter. In fact, the discovery of a neutron star with ~ 2M@ can rule out some of soft EOSs IjDemorest et al 



2O10h , It is also suggested that information about neutron star interiors can be obtained via neutron star asteroseismology o r 
via possible detection of the spect r a of gravitational waves emitted from neutron stars (e.g., lAndersson fc Kokkotas! (1996) 



Sotani. Tominaga fc Maeda ( 2001 ): Sotani. Kohri fc Harada 1 20041 ) : Sotani et al. 1 201lh ). Unlike gravitational waves, there 



are some observational evidences for neutron star oscillations, namely, quasiperiodic oscillations (QPOs) in giant flares from 
soft gamma repeaters (SGRs). Since SGRs are considered to be magnetars, which are neutron stars with strong surface 
magnetic fields of order B > 10 14 -10 15 G, the observed QPOs come presumably from neutron star oscillations, irrespective of 
what portion of a star is actually oscillating. So far, three giant flares have been detected from SGR 0526—66, SGR 1900+14, 
and SGR 1806—20, respectively. The timing analysis of the X-ray afterglow of the giant flares shows that in two o f them, 
specific oscillation frequencies arose in the range from tens of hertz up to a few kilohertz ( Watts fc Strohmaver 20061 ). 
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After this discovery, many theoretical attempts to explain the obs ervations have been made in te r ms of shear torsional oscil- 
lation s in the crust of a neutron star and / or magnetic o scillations (e.g., Levinl d2006 . 2007 ) ; Lee J2007I); Samuelsson fc Andersson 
( 2007h ; ISotani, Kokkotas fe Stergiouiasl (|2007ll2008al ); ISotani, Colaiuda fc Kokkotaj j2008bl ); ISotani fc Kokkotas! j2009l )V The 
important conclusion obtained from such attempts is that either the torsional oscillations or magnetic oscillations dom- 
inate the excited oscillations nea r the surface of a magnet ized neutron star, depending on the magnetic fi eld strength 
I Colaiuda fc Kokkotasll201ll . I2OI2I : iGabler et aDboill . l2012al ) and magnetic configuration jGabler et al. l2012bl ). Given the 
magnetic field strength inferred from the spindown observations ( Kouveliotou et al. 19981 ; Hurley et al.l 19991 ). we can rea- 
sonably consider the QPOs observed in giant flares as crustal torsional oscillations, especially in the case of low-lying modes 
i Passamonti fc Lander! 20131 ). This identification would allow one to probe the properties of ma tter in the crust, which are in 
turn related to the poorly known density depen dence of the symmetry energy of nuclear matter ( Steiner fc Wattal2009l ; Sotani 



2011 ; Gearheart et al. 2011 ; Sotani et al. 20121 ). Once the density dependence of the symmetr y energy is fixed, furthermore, 



one c an address the thickness of the region where nu clei of nonuniform pastalike structures ( Lorenz et al. 19931 ; Oyamatsu 
19931 ) occur in a neutron star (jOvamatsu fc I ida 2007 ^ 



The deeper inside a neutron star, the more uncertainties in the EOS of neutron star matter. Between an ionic ocean 
near the star's surface and a fluid core, a crust occurs. In most of the crust, nuclei are believed to form a bcc Coulomb 
lattice in a roughly uniform electron sea. When the density becomes higher than about 4 x 10 11 g cm~ 3 , neutrons start to 
drip out of the nuclei, and some of them are expected to behave as a superfluid as long as the t emperature is below the 
critical temperature. In fact, such superfluidity plays an important role in modeling pulsar glitches (ISauls!ll989h . while most 
of observed neutron stars are considered to be so cool that most of t he dripped neu trons would behave as a superfluid in 
the absence of the lattice. According to the recent band calculations by Chamel 1 20121 ). however, a significant fraction of the 
dripped neutrons can be entrained non-dissipatively by protons in the nuclei via Bragg scattering off the lattice even at zero 
temperature. Meanwhile, neutron superfluidity due to the remaining dripped neutrons affects the crustal torsional oscillations, 
because such oscillations would be controlled by the enthalpy density of the constituents that comove with the protons 
( van Horn fc EpsteinlHiiol ) . Although there are earlier publications that examine the influence of neutron superfluidity on the 
torsional oscillations in the crust for a specific crust EOS in the Newton ian framework (jAndersson. Glampedakis fc Samuelsson 
20091 : ISamuelsson fc Anderssonll2009l ; IPassamonti fc Andcrss onll2012l). sy stematic calculations in the relativistic framework 



have been recently performed and briefly presented (jSotani et al . 2013a). In this paper, we thus present full details of such 
calculations, which give possible constraints on the parameter L characterizing the density dependence of the symmetry energy 
by comparing the fundamental frequencies of the shear torsional oscillations with the QPO frequencies observed in SGRs. 

The paper is constructed as follows. In the next section, we describe the crust EOSs and equilibrium configurations that 
will be used for calculations of the crustal oscillation frequencies. In the third section, we give the perturbation equation that 
governs the shear torsional oscillations, as well as the boundary conditions, and then address how to take into account the 
effect of neutron superfluidity. In the fourth section, we show the numerical results for the oscillation frequencies and possible 
constraints on L. Finally, the paper closes with a conclusion. We adopt units of c = G = 1, where c and G denote the speed 
of light and the gravitational constant, respectively. 



2 CRUST EQUILIBRIUM CONFIGURATION 

We begin with the properties of matter in neutron star crusts, which are to some degree constrained by empirical data for 
nuclear masses and radii. Near the saturation point of symmetric nuclear matter at zer o temperature, the bulk energy per 
nucleon can be written as a function of baryon number density n b and neutron excess a l|Lattimerlll98ll ): 



w + 



Ko 



(n b - n ) 2 + 



So + - — (n b 

Sn 



- n 0/ 



(1) 



where wo, no, and Ko denote the saturation energy, saturation density, and incompressibility of symmetric nuclear matter. 
Meanwhile, So and L are the parameters associated with the symmetry energy coefficients S(n b ), i.e., So = S(no) and 
L = 3no(dS/dnt) at n b = no- Among these five parameters, Wo, no, and So are fairly well constrained by empirica l masses and 
radii of stable nuclei, while the remaining two parameters L and Ko are left uncertain ( Oyamatsu fc Iida 20031 ). To see how 
strongly the parameters are constrained, two of us (K.O and K.I) first constructed the model for the bulk energy w(n n ,n p ) of 
nuclear matter as a function of neutron and proton number densities, n n = n b (l + q)/2 and n p = nt,(l — q)/2, in such a way 
as to reproduce Eq. ([TJ in the limit of rib no and a — > 0. Then, within a simplified version of the extended Thomas- Fermi 
theory in which the energy density functional for a nucleus of neutron number N and proton number Z, 



E = E h 



■E g +Ec + Nm n + Zm p , 



(2) 



with the bulk energy 



d rn{v)w [n„(r), n p (r)] , 



(3) 
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Table 1. The EOS parameters adopted in the present analysis and the 
corresponding lower and upper densities of the pasta region. 



y (MeV fm 3 ) 


K (MeV) 


L (MeV) 


m (fm 3 ) 


n 2 (fm- 3 ) 


-220 


180 


52.2 


0.060 


0.079 


-220 


230 


73.4 


0.064 


0.073 


-220 


280 


97.5 


0.067 


0.068 


-220 


360 


146.1 


0.066 


0.066 


-350 


180 


31.0 


0.058 


0.091 


-350 


230 


42.6 


0.063 


0.086 


-350 


280 


54.9 


0.067 


0.083 


-350 


360 


76.4 


0.072 


0.076 


-1800 


180 


5.7 


0.058 


0.134 


-1800 


230 


7.6 


0.058 


0.127 


-1800 


360 


12.8 


0.058 


0.118 



the gradient energy 



E„ = F I d 3 r|Vn(r)| 2 , 



the Coulomb energy 



Ec = 



i3 

a r 



3 , np(r)np(r') 

i 1 1 ' 
r — r' 



(4) 



(5) 



and the neutron (proton) rest mass m n {m p 
are for simplicity parametrized as 

1 3 



n;(r) 



1 



(£)" 



o. 



r < Ri, 



r ^ Ri, 



was optimized with respect to the density distributions n n (r) and n p (r), which 



(6) 



where i — n and p. Finally, for given y = —KoSo/(SnoL) and Kq, the most relevant values of wq, no, 5*o, and Fq were obtained 
by fitting the charge number, mass excess, and charge radius, which can be calculated from the optimal density distribution, 
to the empirical valu es for stable nuclei. We remark that the parameter y corresponds to the slope of the saturation line in 
the vicinity of a = |Ovamatsu fc Iidall2003l ). 

In order to ob tain the eq uilibriu m nuclear shape and size as well as the crust EOS for various sets of y and Ko at zero 
temperature, as in Ovamatsu ( 19931 ). Qvamatsu fc Iida (2007) generalized Eqs. ([2]) and to include dripped neutrons of 
uniform number density n^ ut , a neutralizing background of electrons of uniform number density n e , and the lattice energy 
within a Wigner-Seitz approximation. As a result of optimization of the total energy density with respect to the parameters 
characterizing the nucleon distribut ions for given y a n d Kq , the optimal energy density p and nucleon distributions were 
obtained as a function of ni,. As in lOvamatsu fc Iidal (|2007l ); ISotani et al. (2012, l2013al ). we here confine ourselves to the 
parameter range < L < 160 MeV, 180 MeV < Kq ^ 360 MeV, and y < -20 MeV fm 3 , which equa lly well reproduce the 
mass and radius data for stable nuclei and effectively cover even extreme cases jOvamatsu fc Iida| [2003 ) . The EOS parameter 
sets adopted in the present analysis are tabulated in Table [1] where the baryon number densities ni at which the nuclear 
shape changes from sphere to cylinder and ri2 at which the nuclear matter becomes uniform are also li sted. The interval 
betw een ni and 722 corresponds to the pasta region, which decreases with L and vanishes at L ~ 100 MeV (rOvamatsu fc Iida 



2007). We remark that, in order to fill gaps in the values of L, which appear in Table I in ISotani et al.l (|20ia ). we add two 
more parameter sets, namely, (y,K ) = (-220 MeV fm 3 ,280 MeV) and (-350 MeV fm 3 ,280 MeV), where the corresponding 
values of L are 97.5 and 54.9 MeV. 

Let us now consider the equilibrium neutron star configurations. Since the magnetic energy is much smaller than the 
gravitational binding energy even for magnetars, we can neglect the deformation due to the magnetic pressure. Additionally, 
since the magnetars are relatively slowly rotating, we can also neglect the rotational effect. Hereafter, therefore, we consider 
spherically symmetric neutron stars, whose structure is described by the solutions of the well-known Tolman-Oppenheimer- 
Volkoff (TOV) equations. In this case, the metric can be expressed in terms of the spherical polar coordinates r, 8, and <f> 
as 



. 2 2* i,2 . 2A j 2 . 2 , Q 2 , 2 . 2 flJ i2 

as = — e at + e ar +r do + r sm 6 a<p , 



(7) 



where $ and A are functions of r. We remark that A(r) is associated with the mass function 
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i(r) = / dr 4nr' 2 p(r') 



(8) 



2m(r)/r)~ 1 . 

To solve the TOV equations, one generally uses the zero-temperature EOS, i.e., the pressure p as a function of the energy 
density p. For matter in the crust, we use the same EOS models as described above. Unfortunately, the core EOS is still 
uncertain in the absence of clear understanding of the constituents and their interactions both in vacuum and in medium. 
Since we will focus on shear torsional oscillations that occur in the crust, we can effectively describe such uncertainties in the 
core EOS by solely setting the star's mass M and radius R as free parameters, without using specific models for the core EOS. 
In fact, for various sets of M and R, we systematically construct the equilibrium configuration of the crust by int egrating the 
TOV e quations with the c r ust EO S from the star's surface all the way down to the crust-core boundary as in Iida fc Sato! 
(1993); ISotani et all l|2012l . |2013a). This is a contrast to the usual way of constructing a star by initially giving a value of 
the central mass density and then integrating the TOV equations with a specific model for the EOS from the star's center to 
surface. Hereafter, we will consider 1.4 ^ M/Mq ^ 1.8 and 10 km ^ R ^ 14 km as typical values of M and R. Such choice 
of M and R can duly encapsulate uncertainties of the core EOS. 

Generally, a restoring force for shear torsional oscillations is provided by shear stress, which comes from the elasticity 
of the oscillating body and is characterized by the shear modulus [i. In the case of torsional oscillations in the crust of a 
neutron star, the shear modulus is determined by the lattice energy of the Coulomb crystal that constitutes the crust. Since 
the crystal is generally considered to be of bcc type (an fee lattice might occur in place of bec in the crust as suggested by 
recent Thomas-Fermi calculations ( Okamoto et alJl2012D ). one can use the corresponding shear modulus, which is calculated 
for Ze point charges of number density m as 



/i = 0.1194 x 



ni(Ze) 2 



(9) 



where a = (3Z/47rn e ) 1//3 is the radius of a Wigner-Seitz cell ( Strohmaver et al. 1991). Note that this form ula is derived in the 
limit of zero temper ature from Monte C arlo calculations for the shear modulus averaged over all directions ( Ogata fc Ichimarul 
1990h . As shown in ISotani et all (|2012T l, the shear modulus dep ends strongly on the value of L, which comes mainly from 
the L dependence of the calculated Z ( Oyamatsu fc Iidal 120071 ). It is natural that one should take into account the shear 



modulus in pasta p hases, if present, but hereafter we simply assume /i = for the pasta phases, as in lGearheart et al.l (|201ll ) 



IIP _ . 

Sotani et ajj (|2012l . l2013aT ). This is because the shear modulus in the pasta phases except a phase of spherical bubbles has 
at least one direction in which the system is in variant with respect to tra nslation and hence is expected to be significantly 
smaller than that in a phase of spherical nuclei ( Pethick fc Potekhin 19981 ). Under this assumption, we have only to consider 
the shear torsional oscillations that are excited within a crustal region of spherical nuclei, i.e., for n;, ^ n%. Anyway, the 
constraint on L that will be given below can be c onsidered to be robust, because the pasta region is highly limited given the 
resulting constraint on L (jOvamatsu fc ndll2007h . 

We conclude this section by mentioning the thickness of a neutron star crust, i.e., a region of rib ^ ri2. In Fig. [T] we 
illustrate the crust thickness AR obtained from our stellar models f or y = —220 MeV fm 3 . For comparison, we plot the 
parametrization based on the EOS of FPS (|Ravenhall fc Pethi ck 1994): 



(10) 



In this figure, for each value of L, we adopt the 9 stellar models with M = 1.4, 1.6, and 1.8M@ and R = 10, 12, and 14 km. 
From this figure, one can see that the crust thickness in our stellar models is more or less consistent with Eq. (|10p . Although 
we show only the case of y — —220 MeV fm 3 in this figure, the other cases of y — —350 and —1800 MeV fm 3 lap over the 
results for y — —220 MeV fm 3 in Fig.[T] Then, from the 99 stellar models (11 EOS parameter sets in Table[T] 3 different radii, 
and 3 different masses), we derive the fitting formula for the crust thickness as a function of the stellar compactness M/R: 



^„,™(M ) + , 10l( | 



'My 



(11) 



which is also shown in Fig. [T] with thick solid line. We remark that the crust thickness predicted by our stellar models is in 
the range 0.35 km ^ AR ^ 1.35 km. 



3 TORSIONAL OSCILLATIONS 

We now consider the shear torsional oscillations on the equilibrium configuration of the crust of a neutron star described above. 
In order to determine the frequencies, we adopt the relativistic Cowling approximation, i.e., we neglect the metric perturbations 
on Eq. ((7} by setting <5(/ M „ = 0. In fact, one can consider the shear torsional oscillations with satisfactory accuracy even with the 
relativistic Cowling approximation, because the shear torsional oscillations on a spherically symmetric star are incompressible 
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Figure 1. (Color online) Crust thickness calculated for y = —220 MeV fm 3 as a function of the stellar compactness M/R. For 
comparison, the typical behavior suggested by Ravenhall & Pethick (1994) is plotted by the broken line, while the fitting formula l lllH 
is plotted by the thick solid line. 



and thus independent of the density variation during such oscillations. Additionally, due to the spherically symmetric nature 
of the background, we have only to consider the axisymmetric oscillations. Then, the only non-zero perturbed matter quantity 
is the <j> component of the perturbed four- velocity, Su^ , which can be written as 

Su* = e-*d t y(t,r)-±-QdeP e (cose), (12) 

where d t and de denote the partial derivatives with respect to t and 6, respectively, while Pi(cos9) is the ^-th order Legendre 
polynomial. We remark that y (t, r) characterizes the radial dependence of the angular displacement of a matter element. 
By assuming that the perturbation variable y(t,r) has such a harmonic time dependence as y(t,r) = e lut y(r), the per- 
turbation equation that gov erns the shear torsional oscillations can be derived from the linearized equation of motion as 
( Schumaker fc Thorne|[l983r ) 



y" + 



~( 4 ~ + $'-A') + ^ 


y' + 


\r ) fx 





H 2 - 
— uj e 



2A,, 

e y 



0. 



(13) 



where H is the enthalpy density defined as H = p + p with the energy density p and pressure p as described in Sec. [2] and 
the prime denotes the derivative with respect to r. 

Once appropriate boundary conditions are imposed, the problem to solve reduces to an eigenvalue problem with respect to 
uj. Since there is no matter outside the star, we adopt the zero-torque condition at the star's surface. Meanwhile, since there is 
no traction force in the region with p — 0, we adopt the zero-traction condition at the position whe re spherical nuclei disappear 
in the deepest region of the crust. In p ractice, one can show that both conditions reduce to y' = ( Schumaker fc Thorne Il983l ; 
Sotani. Kokkotas fc StergioulasH2007l ). Thus, in determining the frequencies of the shear torsional oscillations, we impose the 
condition of y' = both at n — and ni- 

Now, we take into account the effect of neutron superfluidity on the shear torsional oscillations. In general, it is considered 
that neutrons confined in the nuclei start to drip therefrom when the mass density becomes more than ~ 4 x 10 11 g cm -3 . 
Then, some of the dripped neutrons can behave as a superfluid. Although the behavior of the dripped neutrons is not fully 
understood, a significant fraction of the dripped neutrons may move non-dissipatively with protons in the nuclei as a result of 
Br agg scattering o ff the bcc lattice of the nuclei. In fact, the recent band calculations beyond the Wigner-Seitz approximation 
bv lChameJ (|2012T ) show that the superfluid density, which is defined here as the density of neutrons unlocked to the motion 
of protons in the nuclei, depends sensitively on the baryon density above neutron drip and that a considerable portion of 
the dripped neutrons can be locked to the motion of protons in the nuclei. On the other hand, since the shear torsional 
oscillations are tran sverse, the remaining superfluid n eutrons, whose low-lying excitations are longitudinal, do not contribute 
to such oscillations 1 Pethick. Chamel fe Reddv 2010h . 

We build the effect of neutron superfluidity into the effective enthalpy density H , which can be determined by subtracting 
the superfluid mass density from the total enthal py density H in Eq. (|13p that fully contains the contributions of the superfluid 
neutrons as well as the nuclei and companions ( Iida fc Bavm 20021 ). Since we assume that the temperature of neutron star 



matte r is zero, the baryon chemical potential pb can be expressed as pt — H/n^. Thus, one can write down (jSotani et al 
2013ah 

N, 



H = 



(14) 



where N s denotes the number of neutrons in a Wigner-Seitz cell that do not comove with protons in the nucleus, while A is 
the total nucleon number in the Wigner-Seitz cell. Finally, substituting H for H in Eq. (|13p , one can determine the frequencies 
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Figure 2. (Color online) Effective enthalpy density H calculated as a function of mass density p for y = —220 MeV fm 3 . The solid lines 
are the results for N s /Nd = 0, while the broken lines are for N s /Nd = 1. The four broken lines from top to bottom correspond to the 
cases of L = 52.2, 73.4, 97.5 and 146.1 MeV. 




Figure 3. (Color online) Same as Fig. [2] but for the shear velocity v s . 



of the shear torsional oscillations, which include the effect of neutron superfluidity in a manner that depends on the value of 
N a . Hereafter, we will assume that N 3 comes entirely from the drippe d neutron ga s. Even so, it is still uncertain how much 
fraction of the dripped neutrons behave as a superfluid. Thus, as in ISotani et al.l (|2013al ). we introduce a new parameter 
N s /Nd, where Nd is the number of the dripped neutrons in the Wigner-Seitz cell. For N a /Nd = 0, all the dripped neutrons 
behave as normal matter and contribute to the shear motion, while for N s /Nd = 1, all the dripped neutrons behave as a 
superfluid. We remark that Nd — N s denotes the numb er of the dripp ed neutrons bound to the nucleus. Typically, the value 
of N s /Nd depends on the density inside a neutron star ( Chamelll2012i ). but the case of N s /Nd — in the whole crust is closer 
to the typical behavior than the case of N s /Nd = 1. 

In Fig. [2] we show the effective enthalpies for N s /Nd — by the solid lines and for N s /Nd — 1 by the broken lines, 
where we adopt the EOSs with y = -220 MeV fm 3 and with L = 52.2, 73.4, 97.5 and 146.1 MeV. From this figure, one can 
see that the effective enthalpies for N a /Nd = are almost independent of the EOS parameters, while those for N s /Nd = 1 
depend strongly on the EOS parameters especially for p> 10 13 g/cm 3 . This is because for larger L, the symmetry energy at 
subnuclear densities becomes smaller, leading to increase in the density of the dripped neutrons. Another important quantity 
that characterizes the shear torsional oscillations is the shear velocity v s , which is defined as vl = fi/H. In Fig. [3] we depict 
the shear velocity with the same sets of the EOS parameters, where the solid lines are for N s /Nd = and the broken lines 
are for N s /Nd ~ 1. From this figure, we find that, depending on N 3 /Nd, the shear velocity can double that for N 3 /Nd = 0. 
Both Figs. [5] and |3 suggest the necessity of introducing the effect of neutron superfluidity. 



4 CONSTRAINTS ON THE EOS PARAMETERS 

The shear torsional oscillations are often referred to as t-modes, which are labelled as n te, where £ is the angular index and 
n is the number of radial nodes in the eigenfunctions of the overtones for a specific I. In order to see the dependence of the 
shear torsional oscillations on the EOS parameters, we firs consider the case in which the effect of neutron superfluidity is 
ignored, i.e., N s /Nd = 0. We calculate the fundamental frequencies of such oscillations for a typical neutron star model with 
M = IAMq and R = 12 km by using the 11 EOS parameter sets shown in Table [1] The calculated fundamental frequencies 
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Figure 4. (Color online) The 1 = 2 fundamental frequencies of the shear torsional oscillations, ot2, which are plotted as a function of 
L for M = IAMq and R = 12 km. The thick solid line denotes the fitting formula H15I I. 



(c) 

Table 2. The calculated frequencies, qz^ , for the stellar model of M = IAMq and R = 12 
km in the absence of the effect of neutron superfluidity, and the expected values from Eq. H15I I, 



ot^' • The relative errors, which are determined by [otn^ — ottfi)/ot<>*\ 



are also tabulated. 



y (MeV fm 3 ) 


K (MeV) 


L (MeV) 


o4 C) ( Hz ) 


o4 G) (Hz) 


relative error (%) 


-220 


180 


52.2 


21.44 


20.73 


3.31 


-220 


230 


73.4 


25.56 


25.78 


-0.86 


-220 


280 


97.5 


33.07 


33.04 


0.10 


-220 


360 


146.1 


17.15 


16.60 


3.17 


-350 


180 


31.0 


22.63 


22.91 


-1.24 


-350 


230 


42.6 


32.62 


32.43 


0.57 


-350 


280 


54.9 


13.41 


13.02 


2.96 


-350 


360 


76.4 


19.68 


20.16 


-2.43 


-1800 


180 


5.7 


9.29 


9.40 


-1.13 


-1800 


230 


7.6 


15.32 


16.09 


-5.00 


-1800 


360 


12.8 


30.82 


30.86 


-0.14 



with 1 = 2 are shown in Fig. [4] as a function of L. From this figure, as in ISotani et al.l (|2012f l. one can see that the 1 = 2 
fundamental frequency of the shear torsional oscillations is almost independent of the incompressibility Ko, once the stellar 
model is fixed at M = IAMq and R = 12 km. So is it for different stellar models with IAMq < M < 1.8M and 10 km 
^ R ^ 14 km. Thus, we can focus on the L dependence of the 1 = 2 fundamental frequencies of the shear torsional oscillations. 

Since the number of the EOS parameter sets is limited to eleven, we want to see the L dependence in a continuous 
manner. To this end, we derive a fitting formula for oi2 by assuming the polynomial function form 



oi2 



„(°) 



Co L + Co L , 



(15) 



where Cj ' , c^ 1 ' , and c^ are the adjustable positive parameters that depend on M and R. In practice, we adopt the Levenberg- 
Marquardt algorithm to derive the coefficients Cj ' , c^ 1 ' , and Cj 2 ' . The obtained fitting formula is also shown in Fig. [4] with 
thick solid line. Additionally, for the stellar model with M = 1.4M0 and R = 12 km, we list the calculated 1 = 2 fundamental 
frequencies, ot^ , the expected values from the fitting formula (|15l) . ot^, and the relative errors defined as (o^' — ot^)/ot^ 
in Table [2] For other stellar models, the relative errors are similar to the case of M = 1.4Mq and R = 12 km. This means 
that the fitting formula is in good agreement with the calculated frequencies at least within the accuracy of ~ 5%. 

In Fig. [5] we illustrate the values of ot2 given by the fitting formula (|15[) for the stellar models with 10 km ^ R ^ 14 km and 
IAMq ^ M ^ 1.8Mq, together with the lowest QPO frequency observed from SGR 1806—20 shown in horizontal dot-dashed 
line. We remark that ot2 decreases with R and M. In Fig. [5] therefore, the upper (lower) boundary of the painted region 
corresponds to the 1 = 2 fundamental frequency for M = IAMq (I.SMq) and R = 10 (14) km. Now, on the assumption 
that the QPOs observed in SGR giant flares come from the crustal torsional oscillations, o*2 should become equal to or 
even lower than the lowest frequency in the observed QPOs, because o^2 is theoretically the lowest frequency among many 
eigenfrequencies of the torsional oscillations. Then, from Fig. [5] we can constrain L as L ^ 47.4 MeV if the central objects 
of SGRs are neutron stars with R ^ 14 km and M ^ 1.8Mq. If the oscillating neutron star is a typical one with R = 10 km 
and M = IAMq, we could make a severer constraint on L as L ^ 76.2 MeV. It should be noticed that we omit the effect of 
the pasta phases on the shear modulus in this analysis. Given the resulting constraint on L, the pasta region is too limited 
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L [MeV] 

Figure 5. (Color online) Expected values of 0*2 from fitting formula (1 1 5 1) for the stellar models with 10 km ^ R ^ 14 km and 
1.4M Q sg M ^ 1.8M S . The horizontal dot-dashed line denotes the lowest QPO frequency observed from SGR, 1806-20 (Watts & 
Strohmayer 2006), while the vertical solid and broken lines correspond L = 47.4 MeV and L = 76.2 MeV. 



to have significant consequence to 0^2- In practice, even if we allow for the effect of the pasta phases, the frequency would 
increase because the shear modulus would effectively increase, i.e., the painted region in Fig. [5] would shift to right. Thus, 
our constraint on L is still satisfied. Moreover, even if we take into account the effect of neutron superfluidity, the obtained 



constr aint on L would hold, because the frequencies would increase due to the effect of neutron superfluidity (jSotani et al 



2013al ). 



Next, we consider the effect of neutron superfluidity on the torsional oscillations and try to constrain L by fitting the 
predicted fundamental frequencies of the shear torsional oscillations with different values of I to the QPO frequencies observed 
in SGRs. To take into account the effect of neutron superfluidity, as mentioned before, one needs to know how much fraction of 
the dripped neutrons behave as a superfluid. In fact, earlier calculations of such a fraction are very limited, while the behavior 
of unbound ne utrons inside a nucleus remains to be fully understood. In this paper we adopt the result for N 3 /Nd derived by 
Chamel (2012) as njj/n£, which is based on the band calculations. According to his data, the value of N s /Nd depends on the 



baryon density and becomes around 10 — 30 % at rib ~ 0.01 — 0.4no. Using such a value of N s /Nd, we have calculated the 
fundamental frequencies of the shear torsional oscillations with different values of I for the stellar models constructed with 
the EOS parameter sets as shown in Table [T] Then, we again find that, for each stellar model, the calculated frequencies show 
negligible dependence on Ko but are sensitive to L as in the case of the calculations of 0^2 ignoring the superfluid effect. 
To express the calculated ote as a continuous function of L, we use the same form of fitting formula as Eq. (|15p : 

o** = 4 0) -c^L + cfi 2 , (16) 

where cji , c^ 1 ', and c^ 2 ' are the adjustable positive parameters that depend on M and R for a specific index £. We find that, 
just like the case of Eq. ()15|) . Eq. (|16|) reproduces the calculated frequencies with sufficient accuracy to be used below for 
reasonable fitting to the observed QPO frequencies. Hereafter, we will thus refer to the expected frequencies from Eq. (|16[) as 
the calculated frequencies for simplicity. For comparison of the calculated frequencies with the observed QPO frequencies, we 
particularly focus on the obs erved QPO frequencies lowe r than 100 Hz, i.e., 18, 26, 30, and 92.5 Hz in SGR 1806—20 and 28, 54, 
and 84 Hz in SGR 1900+14 j Watts fc Strohmaverlbood ). because the higher observed frequencies would be easier to explain in 
terms not only of multipolar fundamental and overtone frequencies of shear torsional oscillations in the crust, but also of polar 
type oscillations. Additionally, the possibility to explain the higher QPO frequencies in terms of shear torsional oscillations in 
the had ron-quark mixed phase that may occur in the core of a neutron star is also suggested bv lSotani. Maruvama fe Tatsumi 
1 2013b! ) . 

Due to the small interval between the observed frequencies 26 and 30 Hz in SGR 1806—2 0, theoretical explanations of the 
QPO frequencies observed in SGR 1806—20 are more difficult than those in SGR 1900+14 (jSotani. Kokkotas fc Stergioulas 



2007). Therefore, we first try to reproduce the QPOs observed in SGR 1806—20 by the crustal shear modes. To find a full 



correspondence of the shear torsional oscillations to the obs erved QPOs, o ne should ide n tify th e lowest frequencies in SGR 
1806—20 (18 Hz) as the I = 3 fundamental frequency as in ISotanil (|201ll ); ISotani et al.l (|2013al ). Then, one can manage to 
explain 26, 30, and 92.5 Hz in terms of the fundamental frequencies with I = 4, 5, and 15. In Fig. [6] we compare the predicted 
frequencies with the QPO frequencies observed in SGR 1806—20 for a typical neutron star model with M = IAMq and 
R — 12 km. The best value of L to reproduce the observed frequencies is L — 128.0 MeV, for which the calculated frequencies 
with the best value of L and their relative errors from the observed values are shown in Table [3] We note that the relative 
errors are of the order of such relative errors involved in using the fitting formula (|16[) as shown in Table [2] After performing 
a similar analysis for different stellar models, we find that the QPO frequencies observed in SGR 1806—20 can be explained 
in terms of the eigenfrequencies with the same multipole fundamental oscillations even for the different stellar models, while 
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Figure 6. (Color online) Comparison of the calculated fundamental frequencies of the shear torsional oscillations (solid lines) with the 
QPO frequencies observed in SCR 1806—20 (dot-dashed lines), where we adopt the stellar model with M = 1.4Mq and R = 12 km and 
the Chamel data for Ns/N^. The vertical line denotes the value of L that is consistent with the observations. 



Table 3. The QPO frequencies observed in SCR 1806-20 
and the calculated frequencies with the best value of L to 
reproduce the observed values for the stellar model with M = 
1.4M Q and R = 12 km. 



QPO frequency (Hz) i gt£ (Hz) relative error (%) 

18 3 18.50 -2.79 

26 4 24.82 4.53 

30 5 30.96 -3.19 

92.5 15 90.18 2.51 



the corresponding relative errors are similar to those shown in Table [3] For each stellar model, the obtained best value of L 
is shown in Fig. [7] Assuming that the mass and radius of the oscillating neutron star are in the range 1.4 ^ M/Mq ^1.8 
and 10 km ^ R ^ 14 km, one can constrain L as 101.1 MeV ^ L ^ 160.0 MeV from the observed QPO frequencies of SGR 
1806-20. 

On the other hand, the low-lying QPOs observed in SGR 1900+14 can be similarly explained in terms of the fundamental 
frequencies of the shear torsional oscillations with £ = 4, 8, and 13. For the stellar model with M = 1.4M0 and R — 12 km, 
we show, in Fig. [SI the calculated frequencies as a function of L and compare them with the QPO frequencies observed in 
SGR 1900+14. The best value of L to reproduce the observed frequencies in SGR 1900+14 is L = 113.5 MeV, for which the 
calculated frequencies and relative errors are shown in Table [4] We can likewise obtain the best value of L to reproduce the 
observed QPO frequencies in SGR 1900+14 for each stellar model, which is plotted in Fig. [9] Assuming again that the mass 
and radius of the oscillating neutron star are in the range 1.4 ^ M/Mq ^ 1.8 and 10 km ^ R ^ 14 km, one can constrain L 
as 90.5 MeV a< 131.0 MeV from the observed QPO frequencies of SGR 1900+14. 




Figure 7. (Color online) Values of L at which the calculated fundamental frequencies of the shear torsional oscillations agree best with 
the QPO frequencies observed in SGR 1806—20. The circles, diamonds, and squares correspond to the stellar models with R = 10, 12, 
and 14 km, respectively. 
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Figure 8. (Color online) Same as Fig. [6] but for comparison with the QPO frequencies observed in SGR 1900+14. 



Table 4. Same as Table [3] but for the QPO frequencies 
observed in SGR 1900+14. 



QPO frequency (Hz) 


I 


oh (Hz) 


relative error (%) 


28 


4 


27.26 


2.63 


54 


8 


53.76 


4.50 


84 


13 


86.18 


-2.60 



Note that both observations in SGR 1806—20 and in SGR 1900+14 should simultaneously be explained with a common 
value of L. We can thus obtain a more stringent constraint on L as 101.1 MeV ^ £ ^ 131.0 MeV, which is shown in Fig. 1101 
Some additional remarks are in order. First, such constraint on L in turn would constrain the masses and radii of the central 
objects in SGR 1806—20 and in SGR 1900+14 in such a way that the stellar models out of the painted region are ruled out. 
Such mass-radius constraints would be of particular use in the absence of empiri c al inf ormation about magnetar masses and 
radii. Second, within the nuclear model used here, as shown in Oyamatsu &: Iida 1 20031 ). the symmetry energy at density no, 
So, can be approximately written as a function of L: 



S = 28 MeV + 0.075L. 



(17) 



With this relationship, the constraint on L would translate into 35.6 MeV So ^ 37.8 MeV. One might thus be able to make 
severer constraints on L and So with the help of future possible measurements of magnetar masses and/or radii. 

Meanwhile, as an alternative p ossibility, one ca n identify the three QPOs observed in SGR 1900+14 as the I = 3, 6, 
and 9 fundamental oscillations as in Steiner fc Wattsl |2009l) . As shown in Fig. [TT] this identification is reasonable and indeed 
seems to be in better agreement with the observations than the former one as the I = 4, 8, and 13 oscillations (see Table 0, 
although improvement in the degree of agreement, which is small compared with the intrinsic errors of the fitting formula (|16|l 
used here, may not have to be taken seriously. It is nonetheless interesting to note that the optimal values of L determined by 
the identification as the I = 3, 6, and 9 oscillations are in the range 58.0 MeV ^ L ^ 89.5 MeV, as long as one assumes the 
stellar models with 1.4 ^ M/Mq ^ 1.8 and 10 km ^ R ^ 14 km. This is because such values of L are too small to become 
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Figure 9. (Color online) Same as Fig. [7] but for the best agreement with the QPO frequencies observed in SGR 1900+14. 
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Figure 10. (Color online) The values of L (painted region) that are allowed simultaneously by both of the QPO observations in SGR 
1806—20 and SGR 1900+14 on the assumption that the corresponding neutron stars have mass and radius in the range 1.4 ^ M/Mq ^1.8 
and 10 km ^ R ^ 14 km. 




Figure 11. (Color online) Alternative possible correspondence between the calculated fundamental frequencies of the shear torsional 
oscillations and the QPO frequencies observed in SGR 1900+14. The vertical line denotes the value of L that is consistent with the 
observations. 



consistent with the optimal values of L based on the identification of the four QPOs in SGR 1806—20 as the I = 3, 4, 5, and 
15 oscillations. 

To avoid su ch inconsistency, o ne co uld identify the 18, 30, and 92.5 Hz QPOs in SGR 1806—20 as the I = 2, 3, and 10 
oscillations as in lsteiner fc Watts! (|2009h . although the 26 Hz QPO in SGR 1806-20 remains to be identified. In Fig. [13 we 
show how such an identification works for a typical neutron star model with M = 1.4Mq and R — 12 km; the corresponding 
relative errors are listed in Table HJ] By identifying the low-lying QPOs except the 26 Hz QPO in SGR 1806—20 as the £ — 2, 
3, and 10 shear torsional oscillations, one can obtain the optimal values of L as 54.0 MeV L ^ 85.3 MeV for the stellar 
models with 1.4 ^ M/Mq ^ 1.8 and 10 km ^ R ^ 14 km, which, as in Fig. 1131 ensures the presence of the allowed region 
of L that simultaneously explain both observations in SGR 1806-20 and in SGR 1900+14, i.e., 58.0 MeV < L < 85.3 MeV. 
This region corresponds to 32.4 MeV ^ So ^ 34.4 MeV via Eq. (|17|l . We emphasize that the present identification is again 
in better agreement with the observed frequencies except 26 Hz than the former one as the I = 3, 4, 5, and 15 oscillations. 
This may invoke the possibility that the physics underlying the 26 Hz QPO is missing. We remark that various experim ental 
constraints on L seemingly favor smaller L, although they have yet to converge (see, e.g., Fig. 1 in Newton et al. |2012h V 



Table 5. Same as Table [4] but for the alternative identifi- 
cation. 



QPO frequency (Hz) 


t 


ott (Hz) 


relative error (%) 


28 


3 


27.74 


0.93 


54 


6 


55.48 


-2.74 


84 


9 


82.29 


2.04 
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Figure 12. (Color online) Alternative possible correspondence between the calculated fundamental frequencies of the shear torsional 
oscillations and the QPO frequencies observed in SGR 1806—20 except 26 Hz. The vertical line denotes the value of L that is consistent 
with the observations. 



Table 6. Same as Table [3] but for the alternative identifi- 
cation. 

QPO frequency (Hz) i §tg (Hz) relative error (%) 

18 2 18.23 -1.27 
26 

30 3 28.82 3.93 

92.5 10 94.70 -2.38 



5 CONCLUSION 



In this paper, we have investigated the fundamental frequencies of shear torsional oscillations in neutron star crusts, which 
contain superfluid neutrons, for nine sets of the star's mass and radius as well as eleven models of the crust EOS, and 
considered possible relevance of the crustal shear modes to the QPOs observed from only a few SGRs in the afterglow of giant 
flares. We have succeeded in identifying the low-lying QPOs as the shear modes of different £ in two ways (one of which has 
an obvious caveat), leading to two separate allowed regions of the EOS parameter L characterizing the density dependence of 
the symmetry energy via reasonable fitting of the calculated mode frequen cies to the observed QPO frequencies. The present 
results are basically the same as those obtained in our earlier publications ( Sotani et al. 20121 . 2013a ). but the present work is 
more systematic and quantitatively finer. We hope that future work in this direction, together with other empirical constraints 
on L, will eventually reduce the two allowed regions to one. 

In order to obtain a better constraint on L asteroseismologically, however, many questions remain. Among various 
properties of crustal matter that control the shear torsional oscillations, the shear modulus and the superfluid density play an 
especially significant role. Estimates of the shear modulus in the Coulomb crystals of electron-screened, finite-size charges and 
in the liquid-crystalline pasta phases are thus desired, as well as of the local superfluid density througho ut a Wigner-Seitz cell. 
For the bcc lattice of point charges, for example, the electron screening acts to reduce the shear modulus (|Kobvakov fc Pethick 
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Figure 13. (Color online) Same as Fig. 1101 but for the alternative identification. 
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20131 ). leading to decrease in the eigenfrequencies of the shear torsional oscillations typically by 5 %. Meanwhile, shell and 



pairing effects on the nuclear charge, which are ignored in the present work, are exp ected to have some consequence to the 
shear modulus ( Grill. Margueron. fc Sandulescu 2011 ; Deibel. Steiner. fc Brown 20131 ) . 
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